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ABSTRACT 

We follow the galaxy stellar mass assembly by morphological and spectral type in the COS- 
MOS 2-deg 2 field. We derive the stellar mass functions and stellar mass densities from z = 2 to 
z = 0.2 using 196,000 galaxies selected at F^^^m > lfiJy with accurate photometric rcdshifts 
( <T (zphot-z S pec)/(i+z sp ec) = 0-008 at i + < 22.5). Using a spectral classification, we find that z ~ 1 
is an epoch of transition in the stellar mass assembly of quiescent galaxies. Their stellar mass 
density increases by 1.1 dex between z = 1.5 — 2 and z = 0.8 — 1 (At ~ 2.5 Gyr), but only by 
0.3 dex between z = 0.8 — 1 and z ~ 0.1 (At ~ 6 Gyr). Then, we add the morphological infor- 
mation and find that 80-90% of the massive quiescent galaxies (log(Ai) ~ 11) have an elliptical 
morphology at z < 0.8. Therefore, a dominant mechanism links the shutdown of star formation 
and the acquisition of an elliptical morphology in massive galaxies. Still, a significant fraction of 
quiescent galaxies present a Spi/Irr morphology at low mass (40-60% at log(M) ~ 9.5), but this 
fraction is smaller than predicted by semi-analytical models using a "halo quenching" recipe. We 
also analyze the evolution of star-forming galaxies and split them into "intermediate activity" 
and "high activity" galaxies. We find that the most massive "high activity" galaxies end their 
high star formation rate phase first. Finally, the space density of massive star- forming galaxies 
becomes lower than the space density of massive elliptical galaxies at z < 1. As a consequence, 
the rate of "wet mergers" involved in the formation of the most massive ellipticals must decline 
very rapidly at z < 1, which could explain the observed slow down in the assembly of these 
quiescent and massive sources. 

Subject headings: galaxies: luminosity function, mass function — galaxies: evolution — galaxies: forma- 
tion 



1 Based on observations with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by AURA Inc, under NASA 
contract NAS 5-26555. Also based on observations made 
with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Technol- 
ogy, under NASA contract 1407. Also based on data col- 
lected at : the Subaru Telescope, which is operated by the 
National Astronomical Observatory of Japan; the XMM- 
Newton, an ESA science mission with instruments and 
contributions directly funded by ESA Member States and 
NASA; the European Southern Observatory under Large 



Program 175.A-0839, Chile; Kitt Peak National Obser- 
vatory, Cerro Tololo Inter-American Observatory and the 
National Optical Astronomy Observatory, which are oper- 
ated by the Association of Universities for Research in As- 
tronomy, Inc. (AURA) under cooperative agreement with 
the National Science Foundation; and the Canada-France- 
Hawaii Telescope with MegaPrime/MegaCam operated as 
a joint project by the CFHT Corporation, CEA/DAPNIA, 
the NRC and CADC of Canada, the CNRS of France, TER- 
APIX and the Univ. of Hawaii. 

2 Institute for Astronomy, 2680 Woodlawn Dr., Univer- 
sity of Hawaii, Honolulu, Hawaii, 96822 

3 Laboratoire d'Astrophysique de Marseille, Universite 



1 



1. Introduction 

A clear and comprehensive picture describing 
the physical processes which regulate stellar mass 
growth in galaxies is still missing in our under- 
standing of galaxy evolution. Indeed, the stel- 
lar mass growth is regulated by a complex in- 
terplay between the radiative cooling of the gas 
(e.g. White 1978), cold accretion (e.g. Keres et al. 

2005) , the spatial redistribution of the gas along 
the hierarchical growth of dark matter halos (e.g. 
Springel et al. 2006) and the feedback from super- 
novae and Active Galaxy Nuclei (e.g. Benson et 
al. 2003, Croton et al. 2006). AGN feedback is a 
central process recently added to galaxy formation 
models in order to suppress excessive cooling of the 
gas in massive halos (e.g. Bower et al. 2006, Cro- 
ton et al. 2006, Menci et al. 2006, Cattaneo et al. 

2006) . Even with the inclusion of AGN feedback, 
semi-analytical models still miss a population of 
massive galaxies at z ~ 2 (McCracken et al. 2009) 
and overproduce the number density of low mass 
galaxies (e.g. Kitzbichler & White 2008, Stringer 
et al. 2008). Therefore, a better description of 
star formation activity is still needed. The stel- 
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lar mass function (MF), as studied in this paper, 
characterizes how star formation activity build the 
stellar mass of each galaxy type. 

Merging between galaxies is another central 
mechanism in stellar mass assembly. However, 
there appears to be little consensus between di- 
rect estimates of the merger rate (e.g. Le Fevre et 
al. 2000, Kartaltepe et al. 2007, Lotz et al. 2008). 
An alternative approach is to study the product 
of major mergers. Indeed, these are expected to 
deplete the low mass end of the MF in favor of 
high-mass galaxies, and to produce galaxies with 
elliptical morphologies (e.g. Toomrc & Toomre 
1972). Therefore, a detailed measurement of the 
MF by galaxy type can yield valuable clues on 
galaxy assembly by mergers. This measurement 
can also be considered to be a crucial test of the 
hierarchical paradigm since the assembly of ellip- 
tical galaxies is expected to follow a hierarchical 
build-up similar to that of their host dark matter 
halos (e.g. Kauffmann et al. 1993, de Lucia et al. 
2006). 

Following the stellar mass assembly of a given 
galaxy population requires that the sample be split 
into well characterized galaxy types. A multi- 
color classification scheme is often the only pos- 
sible method to split the faint high redshift sam- 
ples by type. The bimodal distribution of the 
galaxies in a color - magnitude diagram is a com- 
mon tool often used to differentiate two popu- 
lations: "blue cloud" and "red sequence" galax- 
ies (e.g. Bell et al. 2004, Faber et al. 2007, 
Franzetti et al. 2007). The red sequence galax- 
ies include mostly passive galaxies with an el- 
liptical morphology (e.g. Strateva et al. 2001, 
Cassata et al. 2007), but also a significant frac- 
tion of dust extinctcd star-forming galaxies (e.g. 
William et al. 2008) and Spi/Irr galaxies with a 
quenched star- formation (e.g. Bell et al. 2008). 
A novel color-color selection technique (Mjj — My 
versus My — Mj) has been proposed by William 
et al. (2008). This color-color selection breaks the 
degeneracies between dust-extincted star-forming 
galaxies and those with quenched star-formation. 
This diagram is more efficient for detecting a bi- 
modal distribution than a color-magnitude plot 
(William et al. 2008). An alternative multi-color 
classification method is based on a template-fitting 
procedure (e.g. Lin et al. 1999, Wolf et al. 2003, 
Zucca et al. 2006). The advantage of this method 
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is that it defines more than two spectral types. 
But the different template selections are difficult 
to compare from one study to another. 

However, the spectral classifications are sen- 
sitive to the instantaneous star formation rate 
(SFR). Different galaxy populations mixed in the 
same spectral class can be disentangled by adding 
morphological information. Automatic morpho- 
logical classifications (e.g. Abraham et al. 1996) 
performed on high resolution images are efficient 
for discriminating at least two robust classes: 
E/SO and Spi/Irr galaxies (e.g. Lauger et al. 2005, 
Menanteau et al. 2006, Lotz et al. 2008, Capak 
et al. 2008). The combination of morphological 
and spectral classifications allow us to isolate the 
"blue elliptical" galaxies (e.g. Cross et al. 2004, 
Menanteau et al. 2006, Ilbert et al. 2006a) which 
could include newly-formed ellipticals still har- 
boring star-formation (e.g. Van Dokkum & Franx 
2001), Spi/Irr with quenched star formation (Bell 
et al. 2008) and passive elliptical galaxies (e.g. 
Abraham et al. 2007). 

Stellar mass assembly in galaxies by spectral 
and morphological type has already been investi- 
gated using deep optical and near-infrared (NIR) 
surveys. Bundy et al. (2005), Franceschini et al. 
(2006) and Pannella et al. (2006) have derived 
the MF by morphological type using respectively 
2150, 1478 and 1645 galaxies at z < 1.4 in the 
two GOODS fields covering 160 arcmin 2 each (Gi- 
avalisco et al. 2004). Borch et al. (2006) and 
Bundy et al. (2006) derived the MF for blue cloud 
and red sequence galaxies using larger fields of 
0.8 deg 2 and 1.5 deg 2 , respectively. These an- 
alyzes showed that massive elliptical or red se- 
quence galaxies are already in place at z ~ 1, 
while the density is still increasing at lower masses. 
Vergani et al. confirmed these results using the 
4000A Balmer break to separate galaxy popula- 
tions in early and late type systems. Therefore, 
the "downsizing" pattern found by Cowie et al. 
(1996) could be extended to the assembly process 
of ellipticals at z < 1 (e.g. Cimatti et al. 2006). 
Using the K-band luminosity function rather than 
the MF, Arnouts et al. (2008) and Cirasuolo et 
al. (2007) were able to study the stellar mass as- 
sembly for red sequence and blue cloud galaxies at 
z > 1. They found a rapid rise in the space den- 
sity of massive red sequence galaxies from z ~ 2 
to z ~ 1 (Cirasuolo et al. 2007 and Arnouts et 



al. 2007). Abraham et al. (2007) combined mor- 
phology and colors to study stellar mass evolu- 
tion for 144 galaxies with spectroscopic redshifts 
at z > 0.8 using data from the Gemini Deep 
Deep Survey. They confirmed the importance of 
this redshift range in the birth of passive elliptical 
galaxies. 

This paper presents the evolution of the galaxy 
stellar mass function and stellar mass density us- 
ing the COSMOS survey. This survey (Scoville 
et al. 2007) provides four main advantages over 
previous studies that have attempted to measure 
MF evolution: 1) it covers 2-deg 2 which reduces 
the effect of cosmic variance; 2) a morphological 
classification can be carried out based on the Hub- 
ble Space Telescope- Advanced Camera for Surveys 
(HST/ ACS) images (Koekemoer et al. 2007); 3) 
deep Spitzer/IRAC (3.6 — 8.0^m) (Sanders et al. 
2007) and CFHT / WIRC AM i^-band data (Mc- 
Cracken et al. 2009) allow us to estimate accu- 
rate stellar masses out to z ~ 2; 4) the exten- 
sive multi-A coverage of COSMOS provides accu- 
rate photometric redshifts (Ilbert et al. 2009) that 
can be used to derive the galaxy stellar MF. We 
took special care to characterize the galaxy pop- 
ulations, including galaxy morphologies. A first 
study by Scarlata et al. (2007) in the COSMOS 
field already combined morphological and spec- 
tral classifications to study the B-band luminosity 
function. We supplement this study by deriving 
the stellar MF. We provide an estimate of the MF 
which simultaneously covers a large range of red- 
shift (0.2 < z < 2) and a large range of stellar 
masses (10 9 < M/Mq < 10 12 ) using K-band im- 
ages that are 1.5 mag deeper than those used by 
Bundy et al. (2006). We also combine morpho- 
logical and spectral classifications over a field 20 x 
and 100 x larger than Bundy et al. (2005) and 
Abraham et al. (2007), respectively. 

The COSMOS data are introduced in The 
criteria used to split the galaxy sample into various 
populations are described in $3] We present the 
method used to compute the galaxy stellar masses 
in EH $2 SgJ 33 and g8] present the stellar MF 
and stellar mass density of total, early and late 
type samples, respectively. The results are dis- 
cussed in ^21 Throughout this paper, we use the 
standard cosmology (O m = 0.3, Ha — 0.7) 
with Hq = 70 km s _1 Mpc -1 ). Magnitudes are 
given in the AB system. The stellar masses are 
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given in units of solar masses (Mq) for a Chabrier 
initial mass function (hereafter IMF). The stellar 
masses based on a Salpeter IMF (Arnouts et al. 
2007), "diet" Salpeter IMF (Bell et al. 2008) and 
Kroupa IMF (Borch et al. 2006) were converted 
into a Chabrier IMF by adding -0.24 dex, -0.09 
dex, and dex, respectively, to the logarithm of 
the stellar masses. 

2. Data 

2.1. The 3.6/im selected catalogue 

This analysis is based on a mass selected sample 
as generated from the 3.6/im IRAC catalogue of 
the S-COSMOS survey (Sanders et al. 2007). 

The IRAC data were taken during the Spitzer 
Cycle 2 S-COSMOS survey, which used 166 hrs 
to map the full 2-deg 2 COSMOS field (centered 
at J2000 RA =10:00:28.6, Dec = +02:12:21.0). 
The observations were carried out in 4 channels: 
3.6/im, 4.5/im, 5.6/im and 8.0/im. The data were 
initially processed by the Spitzer Science Center 
(SSC). The raw scientific exposures were flux cal- 
ibrated and corrected for well-understood instru- 
mental signatures using a pipeline described by 
Surace et al. (2005). Once the frame- level im- 
ages were prepared, they were projected onto a 
common tangent projection and coadded using the 
SSC MOPEX software^. The images and the cor- 
responding uncertainty maps were generated for 
each of the 4 channels. 

The source catalogue was extracted using the 
SExtractor software (Bertin & Arnouts, 1996). 
The source detection is performed at 3.6/im. The 
IRAC 3.6/im images have a Point Spread Func- 
tion (PSF) of 1.7" which necessitates a careful de- 
blending of the sources. This was obtained with 
a Mexican Hat filtering of the images by SEx- 
tractor. In order to estimate the completeness 
of the 3.6/im catalogue, we simulated point-like 
sources in the 3.6/im mosaic. We simulated si- 
multaneously 10,000 sources with a flux ranging 
from 0.1 to 300 [iJy- The simulated sources were 
distributed randomly in the field without any a 
priory knowledge of the position of the real 3.6/im 
sources (these sources can fall behind or nearby 
a real bright source). We run SExtractor on this 
new image using exactly the same configuration as 
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Fig. 1. — Completeness at 3.6/im: fraction of 
sources simulated in the 3.6/im image which are 
detected with SExtractor, as a function of flux. 
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Fig. 2. — Ratio between the 3.6/im flux measured 
with SExtractor over the simulated flux. 

for real data. Finally, we estimated the fraction of 
simulated sources that we are able to detect, as 
a function of flux (see Figure [T]) . We found that 
the IRAC catalogue is 90% complete at 5 \i Jy and 
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50% complete at 1 \xJy. 

We also used the SExtractor software the mea- 
sure the IRAC fluxes. Following Surace et al. 
(2005), the fluxes were measured over a circular 
aperture of radius 1.9". This small radius pro- 
vides a flux measurement less affected by the pres- 
ence of nearby sources. We tested the accuracy 
of the fluxes recovered by SExtractor using the 
simulation described previously. Figure [2] shows 
the comparison between the simulated flux and 
the flux measured by SExtractor. We obtained 
a flux accuracy of 5%, 10% and 25% for sources 
at 10 < F < 300^Jy, 3 < F < IQ^Jy and 
1 < F < 3fiJy, respectively. These uncertainties 
on the 3.6/irn fluxes are not directly propagated 
into the stellar masses since deep near-infrared 
data (J, H, K) as well as 24 optical bands con- 
strain the rescaling of the best-fit fit templates. 
Based on the same simulation, we derived an aper- 
ture correction of 1.31 at S.Gfim to convert the 
aperture flux to total flux (assuming the sources 
to be point-like). However, we caution the reader 
that specific software like CONVPHOT (De San- 
tis et al. 2006) or TFIT (Laidler et al. 2007) 
could provide flux measurements less affected by 
the confusion by using the K band image as a 
prior. 

Finally, we masked the brightest sources (K s < 
12), as well as poor image quality areas and 
the field boundaries. After removing the masked 
areas, the 3.6^m catalogue contains a total of 
306,000 sources brighter than 1/x Jy (50% complet- 
ness limit) over an area of 2.3-deg 2 . 

2.2. Optical and photo-z catalogues 

We cross-matched the 3.6/ito catalogue with 
the COSMOS photometric (Capak et al. 2008) 
and photo-z catalogs (Ilbert et al. 2009). The 
photo-z were derived for all of the sources in 
the COSMOS photometric catalogue (1,500,515 
sources in total, 937,013 sources at i + < 26.5). 
The photometric fluxes are measured in 31 bands 
(2 bands from GALEX, 6 broad bands from 
SuprimeCam/Subaru camera, 2 broad bands from 
MEGACAM at CFHT, 14 medium and narrow 
bands from SuprimeCam/Subaru, J-band from 
the WFCAM/UKIRT camera, H and K-band 
from the WIRCAM/CFHT camera, and the 4 
IRAC/Spitzer channels). The imaging data are 
extremely deep, reaching u* ~ 27, i + ~ 26.2 and 



K s ~ 23.7 for a 5a detection in a 3" aperture 
(the sensitivities are listed in Capak et al. 2008 
and Salvato et al. 2009). We restricted this study 
to the area covered by the deep optical Subaru 
image (2-deg 2 , 149.4114 < a < 150.8269 and 
1.4987 < S < 2.9127) in order to assure a robust 
photo-z estimate. 

We derived photometric redshifts using the Le 
Phar^\ code (Arnouts et al. 2002 and Ilbert et 
al. 2006b) with a \ 2 template-fitting method. 
The photo-z have been updated in comparison 
to Ilbert et al. (2009) by including new H band 
data. The photo-z are estimated using the median 
of the Probability Distribution Function (PDFz) 
rather than the minimum of the y 2 distribution. 
The photo-z were calibrated with 4,148 spectro- 
scopic redshifts at i\ B < 22.5 from the zCOS- 
MOS survey (Lilly et al. 2009, in prep.). The 
comparison between the photometric and spectro- 
scopic redshifts shows that the fraction of out- 
liers (defined as galaxies with iz v hot — z spec )/{l + 
z S pec) > 0.15) is less than 1% and the accuracy 
is as good as <7(z phot -z Bpec )/(.i+z apec ) = 0.008 at 
i\ B < 22.5. A spectroscopic follow-up of 24/im 
selected sources at z < 1.5 (Kartaltepe et al. 
2009, in prep.) allows us to characterize the 
photo-z accuracy at fainter magnitude. We found 
V(z phot -z spsc )/(i+ Zsp e C ) = 0.011 at 22.5 < %\ B < 24 

aIld Cr (2 ph ot-^pec)/(l+^pao) = - 053 &t 24 < ^AB < 

25 for this infrared selected sample. At z > 1.5, 
we used the zCOSMOS-faint spectroscopic sam- 
ple (Lilly et al. 2009, in prep.) to quantify the 
quality of the photo-z in the magnitude/redshift 
range where the photo-z are expected to have the 
highest uncertainty. These color selected galaxies 
have median apparent magnitude of «^ od ~ 24.1 
and a median redshift of 2.1. At 1.5 < z < 3, 
we obtained an accuracy of 0a 2 /(i+ 2 ) = 0.04 with 
10% of catastrophic failures. However, these vari- 
ous spec-z samples probe only specific populations 
(infrared selected, color selected). Figure 8 of Il- 
bert et al. (2009) shows that the photo-z la error 
derived from the PDFz is well representative of the 
photo-z accuracy. The median la error is 0.02 for 
the full catalogue at F 3 , 6 fj.m > l[i-Jy and 0.08 in 
the redshift range 1.25 < z < 2. We also showed 
in Figure 12 of Ilbert et al. (2009) that the photo- 
z accuracy is degraded at i\ B > 25.5. Therefore, 
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we take special care in limiting the contribution of 
these faint sources in our analysis (see §4.4|) . 
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Fig. 3. — Photometric redshift distributions 
for the 3.6/im selected sample (ttiab (3.6/im) = 
-2.5/og(F)+23.9 where F is the flux in fiJy). The 
dotted line corresponds to the redshift distribution 
without the galaxies not detected in optical. 

2.3. Identification of the 3.6/im source 
counterparts 

We cross-matched the 3.6/xm and photo-z cat- 
alogues by taking the closest counterpart within a 
radius of 1". The distances of the first and second 
closest optical counterpart have a median value of 
0.2" and 2.5", respectively. The two distributions 
intersect at 1". Therefore, we adopted a matching 
distance of 1" which is a good compromise to de- 
tect the maximum of optical counterparts and lim- 
its the risk of wrong identification. Still, 2.6% of 
the IRAC sources have two possible optical coun- 
terparts in less than 1". In order to estimate the 
probability of having identified the wrong opti- 
cal counterpart, we multiplied the probability of 
having the right counterpart by the probability of 
having another optical source at a lower distance 
(using the distance distributions of the first and 
second closest counterpart, respectively). We ob- 
tained that the probability of having identified the 
wrong optical counterpart is 0.1%. We found a 



similar probability of 0.4% using simulations. 

We identified 8507 3.6/im sources without opti- 
cal counterpart (about 4% of the IRAC catalogue) . 
Most of these sources are extremely faint at 3.6/im, 
without counterpart in the K-band selected cata- 
logue (McCracken et al. 2009). This sample in- 
cludes also a significant fraction of fake detections 
created by the residual of the muxbleed correction 
(Surace et al. 2005). Still, we were able to identify 
2714 IRAC sources which are clearly non-detected 
in optical and are detected in the K-band selected 
catalogue. These sources can be 2 > 1.5 quies- 
cent systems. Therefore, we included them in our 
analysis. We measured a photo-z for these sources 
using NIR and IRAC data. An upper-limit was set 
in i + since this band was used for galaxy detection 
in Capak et al. (2009). The averaged redshift of 
this population is z ~ 2.9. 93% and 77% of these 
sources are at z > 1.5 and z > 2, respectively. 
In any case, the impact of this galaxy population 
on our analysis is low since the stellar mass lim- 
its are set to ensure a low fraction of sources with 



i\ B > 25.5 in the stellar mass sample (see g4T4J). 

Finally, we removed all of the sources flagged as 
star or AGN. Stars were removed from the sample 
by comparing the x 2 evaluated for both the galaxy 
templates and stellar templates (see §3.6 of Ilbert 
et al. 2009). The 1,887 sources (1% of the total 
sample) detected with XMM-COSMOS (Hasinger 
et al. 2007, Brusa et al. 2007, Salvato et al. 2009) 
were removed from the sample since their optical 
emissivity is likely dominated by an AGN. 

To summarize, this study is based on the S- 
COSMOS 3.6/xm selected catalogue which is 50% 
complete at 1 fxjy. We cross-match this cata- 
logue with the full optical and photo-z catalogue 
using a match distance of 1". The photo-z accu- 
racy is as good as <T(z phot -z 3peB y(i+z 3peB ) = 0.008 
at i\ B < 22.5. The final sample (after having re- 
moved the stars, the XMM sources, the masked ar- 
eas and the objects without optical counterparts) 
contains 196,000 galaxies at F 3 , 6tim > 1/iJy over 
an effective area of 1.73-deg 2 . Figure [3] shows the 
redshift distributions for the F^.^^m > lfiJy se- 
lected sample with a median redshift of z ~ 1.1. 
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3. Morphological and spectral galaxy clas- 
sifications 

In this paper, we study the galaxy stellar mass 
functions per morphological and spectral type. 
This section presents the criteria used to define 
the various galaxy types: E/SO and Spi/Irr galax- 
ies based on their morphology ; three spectral 
classes ("quiescent", "intermediate activity" and 
"high activity" galaxies) using best-fit templates. 

3.1. Morphological classification 

We used the high resolution HST/ACS images 
(Koekemoer et al. 2007) to perform a morphologi- 
cal classification of our galaxy sample. The images 
in the F814W filter reach a depth of 27.8 mag for a 
point source at 5a. We adopted two independent 
morphological classifications to separate E/SO and 
Spi/Irr galaxies. 

The first classification is based on the Gini (G) 
and concentration (C) parameters measured by 
Abraham et al. (2007) (hereafter G-C classifica- 
tion). The Gini parameter measures the inequal- 
ity with which the light of a galaxy is distributed 
among its constituent pixels. Like Capak et al. 
(2007), the galaxies with G > 0.43 were consid- 
ered E/SO galaxies. In addition, we rejected from 
the E/SO sample the galaxies with a concentration 
parameter smaller than 0.3 (Ilbert et al. 2006a). 

The second classification was performed by 
Cassata et al. (2009, in prep.) (hereafter C09 
classification). The structural parameters are 
measured using a "quasi-Petrosian" image tresh- 
olding technique (Abraham et al. 2007). This 
classification includes Gini, Concentration, Asym- 
metry, and M20 (e.g. Lotz et al. 2004). The 
multi-dimensional parameter space is automati- 
cally converted into an E/SO and Spi/Irr classifi- 
cation by matching these parameters with those of 
a training sample of 250 visually classified galaxies 
(50 galaxies per 0.5 mag bin out to i\ B < 24). 

The E/SO selection performed by C09 is more 
conservative than the G-C classification (i.e. less 
contaminated by spiral galaxies), but is likely to 
be more incomplete. Indeed, less than 1% of the 
E/SO sources from C09 are not identified as E/SO 
with the G-C parameters, while 33% of the E/SO 
sources identified with the G-C parameters are not 
identified by C09 (at log{M) > 10 and z < 1.2). 




2 4 6 



Fig. 4. — Distribution of the unextincted rest- 
frame color (NUV — r + )tem P iate of 1,500 visually 
classified galaxies (late spiral and irregular: blue 
dotted line; early spiral: green dashed line; E/SO: 
red solid line) . The vertical dashed lines show the 
separation between the "quiescent" , "intermediate 
activity" and "high activity" galaxies. 

3.2. Spectral classification 

A set of SED templates was generated using 
the Bruzual & Chariot (2003) (BC03) package 
and fitted to the multi-color data (see section l4~Tjl . 
The extinction is added as a free parameter in 
the fit. We used the unextincted rest-frame col- 
ors (NUV — r + ) template of the templates to define 
three spectral classes: a) the "quiescent" galax- 
ies with (NUV - r + ) t e mp i a te > 3.5; b) the "in- 
termediate activity" galaxies with 1.2 < (NUV — 
r + ) template < 3.5; c) the "high activity" galaxies 
with (NUV - r+) temp iate < 1.2. 

Figure [4] shows the (NUV— r + ) terrap ; ate distribu- 
tion of a sample of 1,500 galaxies that we visually 
classified as E/SO, early spiral, late spiral or irreg- 
ular. The 1,500 galaxies were selected to provide 
an unambiguous visual classification (isolated and 
bright galaxies) but were not selected to be sta- 
tistically representative of the 3.6/xm sample. A 
cut at (NUV — r + )t emp iate > 3.5 isolates well the 
E/SO galaxies. The "intermediate activity" class 
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Fig. 5. — Rest-frame colors M(NUV) — M(r + ) versus M(r + ) — M(J) (not corrected for dust reddening) 
from z = 0.2 (top left panel) to z — 2 (bottom right panel). The red open triangles, green crosses, and blue 
points are the galaxies selected as "quiescent" , "intermediate activity" and "high activity" , respectively, on 
the basis of their unextincted rest-frame color (NUV — r + ) template- 



(1.2 < (NUV - r + ) temp iate < 3.5) includes most 
of the visually selected early spiral galaxies but is 



strongly contaminated by late spiral and irregular 
galaxies. 
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Fig. 6.— Distribution of the Specific SFR (SSFR) 
for the "quiescent" (red dotted line), "intermedi- 
ate activity" (green solid line) and "high activity" 
(blue dashed line) galaxies (at 0.2 < z < 1.2 and 
log(M) > 10) (se 



Figure [5] shows a slightly modified version of the 
color-color selection technique (Mjj — My versus 
My - Mj) proposed by William et al. (2008). We 
used the color NUV — r + instead of U — V since 
this color is a better indicator of the current ver- 
sus past star formation activity (e.g. Martin et 
al. 2007, Arnouts et al. 2007). The rest-frame 
colors were computed as described in appendix [A] 
and are not corrected for internal dust attenuation 
(by contrast with (NUV — r + ) template which is cor- 
rected for dust attenuation) . A red clump appears 
clearly from z = 0.2 out to z = 2. This clump is 
mostly composed of "quiescent" galaxies. There- 
fore, our "quiescent" population is similar to the 
red clump population selected by William et al. 
(2008). The galaxies with a red NUV - r+ > 4.5 
rest frame color are well separated into a "quies- 
cent" population with r + — J < 1.2 (red clump) 
and a dust-extincted star-forming population with 
r+ - J > 1.2. 

Finally, we show in Figure [6] that each spec- 
tral class corresponds to a range of Specific SFR 
(SSFR), computed as the instantaneous SFR from 
the best-fit template divided by the stellar mass. 



To summarize, we used the HST/ACS images 
to separate E/S0 and Spi/Irr galaxies using two 
morphological classification methods (C09 and G- 
C). We also defined three spectral classes on the 
basis of the best-fit templates which are "quies- 
cent" , "intermediate activity" and "high activity" 
galaxies. These three spectral classes are well sep- 
arated in ranges of SSFR. We showed that our 
"quiescent" population matches well with the red 
clump galaxies found by Williams et al. (2008) 
and is consistent with an E/SO population selected 
morphologically. 
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Table 1: Parameters used to generate the SED 
templates with the BC03 package. 



4. The galaxy stellar mass sample 

In this section, we describe the method used 
to measure galaxy stellar masses and the galaxy 
stellar mass function. 

4.1. Technique used for estimating stellar 
masses 

We used stellar population synthesis (SPS) 
models to convert luminosity into stellar mass 
(e.g. Bell et al. 2003, Fontana et al. 2004). 
The stellar mass is the factor needed to rescale 
the best-fit template (normalized at one solar 
mass) for the intrinsic luminosities. The SED 
templates were generated with the stellar popula- 
tion synthesis package developed by (BC03). We 
assumed an universal IMF from Chabrier (2003) 
and an exponentially declining star formation his- 
tory SFR cx e~*/ r (r in the range 0.1 Gyr to 30 
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Fig. 7. — The black solid line histograms show the difference between the stellar masses computed with 
photo-z and spectro-z. We used the zCOSMOS spectroscopic sample selected at i\ B < 22.5 (Lilly et al. 
2009, in prep.) in the left panel and the spectro-z of infrared selected sources from Kartaltepe et al. (2009, 
in prep.) in the right panel. The thin black dashed lines are gaussian distributions with a = 0.02 (left panel) 
and a = 0.03 (right panel). The green dashed lines show the difference between the stellar masses computed 
with BC03 and Chariot & Bruzual (2007, private communication). The magenta dotted lines show the 
differences between the stellar masses computed using the Calzetti et al. (2000) and Chariot & Fall (2001) 
extinction laws. The redshifts were set to the spectro-z values in the two last cases. Systematic uncertainties 
due to the models dominate the errors introduced by the photo-z, at least in the magnitude/redshift range 
explored with our spectroscopic samples. 



Gyr). The SEDs were generated for a grid of 51 
ages (in the range 0.1 Gyr to 14.5 Gyr). Dust 
extinction was applied to the templates using the 
Calzetti et al. (2000) law (E(B-V) in the range 
to 0.5). We used models with two different metal- 
licities. The parameters used to generate the SED 
templates are listed in Table [TJ Following Fontana 
ct al. (2006) and Pozzetti et al. (2007), we im- 
posed the prior E(B - V) < 0.15 if age/r > 4 (a 
significant extinction is only allowed for galaxies 
with a high SFR). 

We introduced the fluxes measured at 24 /urn 
with the Spitzer /MIPS instrument (Aussel et al. 
2009, in prep.) as a constraint in the template 
fitting procedure, as detailed in appendix lEl how- 
ever, the 24 /im constraint had little effect on the 
derived stellar masses. The dispersion between the 
stellar masses with and without using this con- 
straint is 0.014 dex, and only 1% of the sample 



differs by more than 0.2 dex. 

All the available bands were used to compute 
the stellar masses (broad bands as well as medium 
bands). The IRAC data were included in the fit. 
We compared the stellar masses computed with 
and without IRAC. We obtained a dispersion of 
0.02, 0.02, 0.04, 0.11 dex at z=0.2-0.5, z=0.5-l, 
z=l-1.5, z=1.5-2, respectively. At z < 1.5, the 
IRAC data have a small impact on the stellar 
masses, negligible in comparison to the 0.2 dex 
uncertainty expected in the stellar mass estimate 
(e.g. Pozzetti et al. 2007, Longhetti & Saracco 
2008). At z > 1.5, including the IRAC data mod- 
ifies significantly the stellar masses. The IRAC 
data are necessary since the K band doesn't probe 
anymore the NIR rest-frame wavelength range. 
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4.2. Systematic uncertainties in the stellar 
mass estimate 

We quantified how the stellar mass accuracy is 
affected by the use of photo-z rather than spectro- 
z. Figure[7]shows the difference between the stellar 
masses computed with the photometric and spec- 
troscopic redshifts. We used two spectroscopic 
samples: the zCOSMOS bright spectroscopic sam- 
ple selected at i+ B < 22.5 (Lilly et al. 2009, 
in prep.) and a spectroscopic follow-up of 24/xto 
selected sources (median flux F 2 4 Mm ~ lAOfiJy) 
by Kartaltepe et al. (2009, in prep.). The in- 
frared follow-up supplements very well the zCOS- 
MOS spectro-z since the former sample is fainter 
(18 < i\ B < 25 with 43% of the sources be- 
ing fainter than i\ B = 22.5) and extends out to 
z ~ 1.5 (median redshift of z ~ 0.74). We found a 
median difference smaller than 0.002 dex between 
the photo-z and the spectro-z stellar masses for 
both samples. Therefore, no systematic offsets 
appear to be introduced by the use of our photo- 
z. The dispersion between the two estimates is 
smaller than ~ 0.03 dex. This dispersion is 10 x 
smaller than the systematic uncertainties expected 
in the stellar mass estimate (e.g. Pozzetti et al. 
2007, Longhetti & Saracco 2008), at least in the 
magnitude/redshift space covered by the spectro- 
scopic samples. 

The choice of extinction law impacts the fit of 
the template and therefore the mass-to-light ra- 
tio. We computed the stellar masses using both 
the Calzetti et al. (2000) and the Chariot & Fall 
(2000) extinction laws (the latter is included in 
BC03). The redshifts are set to the spectro-z val- 
ues. The median difference between the two stellar 
mass estimates (Calzetti - Chariot & Fall) is -0.14 
dex with a dispersion of 0.10 dex for the zCOS- 
MOS sample (left panel of Figure [7]) . This median 
difference reaches -0.27 dex for the MIPS spectro- 
scopic sample, showing that the systematic offset 
is larger for massive galaxies with a high SFR. We 
adopted the Calzetti et al. (2000) extinction law. 
This choice is favored by a comparison between 
the SFR derived from the best-fit template and 
the SFR measured from mid-infrared 24 fj,m data, 
as described in appendix [B] 

The stellar mass estimate depends also on the 
assumed stellar population synthesis model. We 
computed the stellar masses with an upgraded ver- 



sion of the BC03 model including a better treat- 
ment of the thermally pulsing asymptotic giant 
branch (TP-AGB) phase (Bruzual C, 2007, Char- 
lot & Bruzual, 2007, private communication). Fig- 
ure [7] shows a comparison of the stellar mass es- 
timates using the two versions of the models. We 
find a median difference of 0.13-0.15 dex and a 
dispersion of 0.09 dex. Pozzetti et al. (2007) mea- 
sured a difference of 0.14 dex between the stellar 
masses computed with the BC03 and the Maras- 
ton (2005) models (which also includes treatment 
of the TP-AGB stars) , in agreement with our re- 
sults. We used the public version of BC03 for con- 
sistency with results from the literature. However, 
the MF computed in this paper have also been 
computed with the upgraded version of BC03, and 
our conclusions remained unchanged. 

4.3. Method to estimate the mass function 

We measured the stellar mass functions using 
the tool ALF (Algorithm for Luminosity Func- 
tion, Ilbert et al. 2005). This tool includes the 
STY parametric estimator (Sandage, Tammann 
& Yahil 1979) and three non-parametric estima- 
tors: the l/V ma x (Schmidt 1968), C+ (Lynden- 
Bell 1971 and Zucca et 1997) and the Step- Wise 
Maximum Likelihood (SWML, Efstathiou 1988). 
A brief introduction of these standard estimators 
is given in appendix [Cl 

We weighted each galaxy according to the com- 
pletness of the 3.6/j.m catalogue (see section 12.11 
and Figure [I]). We attributed a weight to each 
source depending on the 3.6/im flux. This weight 
is the inverse fraction of 3.6/im sources detected at 
this flux (e.g. a weight of 2 is given to the sources 
at F 3 . 6fim ~ 1/xJy). 

We also performed extensive simulations in or- 
der to propagate the photo-z uncertainties into 
the mass function. A redshift probability distri- 
bution function (PDFz) was attributed to each 
galaxy when we measured the photo-z (see Ilbert 
et al. 2009). We created 20 catalogues by ran- 
domly picking a redshift within the PDFz of each 
object. The MFs were measured for each of the 
20 catalogues in every redshift bin and for each 
galaxy population. Finally, the dispersion of the 
Schechter parameters was measured over the 20 
realizations. We added in quadrature the Poisso- 
nian errors and these errors induced by the use of 
photo-z. 
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Fig. 8. — Fraction of galaxies with an apparent 
magnitude fainter than i + = 24 (solid lines) and 
i + = 25.5 (dashed lines) as a function of stellar 
mass. This fraction is measured per redshift bin 
(top to bottom panels). The left and right pan- 
els correspond to the quiescent and star-forming 
galaxies, respectively. From this fraction, we de- 
fined the lowest stellar mass limit which ensures a 
maximum of 30% of galaxies fainter than i + > 24 
and i + > 25.5 in the lowest stellar mass bin of the 
MF. 



4.4. Considered stellar mass range for the 
MF estimate 

The low mass limits considered for the MF es- 
timates are set in order to insure a complete and 
unbiased stellar mass sample with accurate photo- 
z. 

We defined the low stellar mass limits in or- 
der to reduce the fraction of optically faint sources 
with low quality photo-z in the stellar mass sam- 
ple. Figure [5] shows the fraction of galaxies with 
> 25.5 (dashed lines) as a function of the stel- 



;+ 



lar mass. We defined the stellar mass ranges in 
order to keep the fraction of galaxies fainter than 
i + > 25.5 below an arbitrary limit of 30%. We set 
the limit to i\ B — 25.5 since the photo-z are de- 
graded at fainter magnitudes (Ilbert et al. 2009). 
With this approach, the lowest stellar mass bin of 



the MF has 30% of its objects with lower accu- 
racy photo-z, and this fraction decreases rapidly 
in the higher stellar mass bins. According to 
Abraham et al. (2007) and Capak et al. (2007), 
the morphological classification is robust out to 
%\ B ~ 24. When a morphological selection is ap- 
plied, we therefore adopted a limit at i\ B < 24 
(solid line in Figured]) rather than i\ B < 25.5. 

Moreover, Ilbert et al. (2004) and Fontana et 
al. (2004) showed that the MF estimators can be 
biased at low masses because galaxies with differ- 
ent SEDs (and mass-to-light ratios) are not visible 
up to the same stellar mass limit. This affects each 
MF estimator differently (see Figure 4 in Ilbert et 
al. 2004). We restricted our MF estimate to the 
stellar mass range where the 3 non-parametric es- 
timators agree, to limit the impact of such bias on 
our results. 

The lowest mass limits considered for the MF 
estimate are given in Table [5] and Table [3] for each 
redshift bin. In all cases, the lowest limit in mass 
has to be the largest value imposed by the mix of 
galaxy type and the photo-z limitation at i\ B < 
25.5. 



quiescent 



intermediate 
activity 



high 
activity 



star-forming 



red elliptical blue elliptical 



elliptical 



spiral/lrr 



Fig. 9. — Schematic view of the different classifi- 
cations used in this paper. The top and bottom 
panels are the spectral and morphological classifi- 
cations, respectively. 
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log(M*) [MJ 

Fig. 10. — Total MF (no separation by type). The open circles are the non-parametric estimates of the MF 
(1/Vmax)- The solid black lines are the sum of the "quiescent", "intermediate activity", and "high activity" 
MFs as taken from Tableland Table [3] The total MFs are compared with data from the literature (all MFs 
are converted into a Chabrier IMF). 



5. Mass Function of the total sample 

We first analyze the total MF (no cut by mor- 
phological or spectral type) and compare it with 
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data from the litterature. In the next sections, we 
will study the MF per morphological and spectral 
type (Figure [9] summarizes the different classifica- 
tions used hereafter.) 

Figure [10] shows the estimate of the total MF. 
The non-parametric MF estimate (open black cir- 
cles) shows a small turn-over at log(A4) < 10 and 
z < 0.8 which can not be reproduced well with a 
Schechter function (Schechter 1976, see appendix 
|C|) . This turn-over is also detected in Pozzetti et 
al. (2009) and Drory et al. (2009). Therefore, 
we do not give a Schechter parametrization of the 
total MF. A parametrization can be retrieved by 
summing the Schechter fits of the "quiescent" , "in- 
termediate activity", and "high-activity" galaxies 
given in Tableland Table[3] This sum (solid black 
curves in Figure [TO]) provides a parametrization in 
excellent agreement with the non-parametric esti- 
mate. 

We compared the total MF and data from the 
literature (all MFs are converted into a Chabrier 
IMF and to the same cosmology). In general, 
we find excellent agreement between the differ- 
ent MFs out to z = 2. The offsets between the 
high-mass exponential cutoffs (i.e. the sharp de- 
cline of the density above the characteristic stellar 
mass M.*) are smaller than 0.2 dex. Combined dif- 
ferences due to cosmic variance and methodology 
used to measure the stellar massed are consistent 
with differences of 0.2 dex (see £|4.2[) . 

6. Stellar Mass assembly of quiescent and 
elliptical galaxies 

In this section, we present the MF for galaxies 
having a "quiescent" spectral type. From the tight 
correlation seen in the local Universe between mor- 
phology and colors, the quiescent galaxies are ex- 
pected to preferentially have an elliptical morphol- 
ogy. However, the correlation between color and 
morphology is not perfect (Bell et al. 2008) and 
we need to quantify how this relation evolves with 
redshift. Therefore, we also derived the MF of 
quiescent galaxies having an elliptical morphology 



3 For instance, Bundy et al. 2006 used the Chariot & Fall 
2001 extinction law, Pozzetti et al. 2007 did not allow sub- 
solar metallicities, and Borch et al. 2006 used the PEGASE 
stellar population synthesis package from Fioc & Rocca- 
Volmerange 1997 and no NIR data. 



( "red elliptical" , see Figure [9]) . 

6.1. MF of quiescent galaxies 

Figure QT] shows the MF of the quiescent galax- 
ies in the range z = 0.2 — 2. We find that: i) 
the density of "quiescent" galaxies more massive 
than log(M.) > 11 increases by a factor of ~ 14 
between z = 1.5 — 2 and z — 0.8 — 1; ii) this evo- 
lution slows down significantly after z < 1 and the 
high-mass exponential cutoff does not increase by 
more than 0.2 dex at z < 1; iii) the density of 
quiescent galaxies increases at intermediate mass 
between z = 0.8 — 1 and z = 0.2—0.4 (e.g. by a fac- 
tor of 4.4 at log(M) ~ 10). In Figure [T2l we have 
over-plotted our results and the local measurement 
performed by Bell et al. (2003). We find consistent 
evolutionary trends when we compare our data to 
the local measurementQ: the local density is higher 
at intermediate masses (9 < log(A4) < 11) and the 
local exponential cutoff is consistent within 0.2 dex 
with the values obtained at Q.2 < z < 1. 

We analyze the evolution of the best-fit Schechter 
parameters of the quiescent MF (Table [2]). The 
top panel of Figure 1131 shows a continuous steep- 
ening of the slope a with time. This steep- 
ening reflects the rapid density increase of the 
low/intermediate mass galaxies. In the middle 
panel, the normalization $* of the quiescent MF is 
shown to increase by a factor of 15 from z — 1.5 — 2 
to z = 0.8 — 1. The rapid increase of <&* is no 
longer detected from z = 0.8 — 1 to z = 0.2 — 0.4 
where remains approximately constant. Some 
fluctuations (a factor of ~ 2) appear in this red- 
shift range when we reduce the size of the redshift 
bins to Az = 0.1, which is consistent with cosmic 
variance (Scoville et al. 2007). Finally, the char- 
acteristic stellar mass M.* increases by 0.3-0.4 dex 
between z ~ 1 and z ~ 0.3 (bottom panel). 

6.2. MF of red elliptical galaxies 

Figure QT] shows the MF of the red ellipticals 
(quiescent with an elliptical morphology) galaxies 
in the range z = 0.2 — 1.2. We find a similar evo- 
lution for the red ellipticals as for the quiescent 
galaxies. Between z = 0.8 — 1 and z = 0.2 — 0.4, 



4 The local MF is computed with the code PEGASE. The 
local MF could be shifted by +0.06 dex to match our stel- 
lar masses computed with the BC03 code (Rettura et al. 
2006). 
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Quiescent galaxies 




log(M.) [M e ] 

Fig. 11. — Evolution of the MFs for "quiescent" galaxies from z = 0.2 — 0.4 (top left panel) to z — 1.5 — 2 
(bottom left panel). The "quiescent" galaxies are shown by black open circles and solid lines. The lower 
and upper envelopes of the shaded area are the red elliptical classified using C09 and G-C, respectively. The 
black long-dashed line is the quiescent MF estimated at z = 0.2 — 0.4; it is shown in each panel to serve as 
a reference. 

the density of red ellipticals increases by a factor high mass (log(M) > 11). Therefore, the most 

of 4-5.4 at intermediate mass (log(M.) ~ 10). By massive red elliptical galaxies show little evolu- 

contrast, their density increases only by 1.7-2 at tion at z < 1 while their density still increases at 
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Fig. 12. — MF of the "quiescent" galaxies from 
z = 1.5 - 2 to z = 0.2 - 0.4. The numbers on the 
MFs are the rcdshifts. The green vertical boxes are 
an estimate of the cosmic variance by Scoville et 
al. (2007) for halo mass ranges of 10 13 - 1O 14 7W . 
The black vertical lines correspond to the two ex- 
treme MFs in 4 sub- fields obtained by splitting the 
2-deg 2 into four quadrants each of 0.5-deg 2 . The 
red solid line and points are the local MF mea- 
surements for red galaxies by Bell et al. (2003) at 
z ~ 0.1. 

low/intermediate masses. 

The ratio between the red elliptical and qui- 
escent MFs is plotted in Figure [TU This ra- 
tio quantifies the fraction of quiescent galaxies 
with an elliptical morphology (and the comple- 
mentary information about the fraction of Spi/Irr 
with quenched star formation). At high mass, 
log(Ai) ~ 11, the fraction of quiescent galaxies 
with an elliptical morphology is greater than 90% 
(80%) at z < 0.8 for a G-C (C09) classification. 
This fraction has a maximum at 10.5 < log(M) < 
11 and decreases continously toward low masses 
reaching 60% (40%) at log(M) ~ 9.5. This frac- 
tion seems also to decrease at really high mass 
(log(A4) ~ 11.5). However, we caution the reader 
that the constraint on the MF at such high masses 
relies on few galaxies. 

The fraction of quiescent galaxies with an el- 



Fig. 13. — Evolution of the Schechter parameters 
as a function of redshift (black solid lines: quies- 
cent; red dashed lines: red ellipticals). The long 
dashed magenta lines show the Schechter parame- 
ters measured in smaller redshift bins of Az = 0.1 
for the quiescent galaxies. From the top to the 
bottom: evolution with redshift of the slope, the 
normalization and the characteristic stellar mass, 
respectively. The vertical lines for the quiescent 
galaxies represent the extreme values in four sub- 
fields of 0.5-deg 2 after having divided the COS- 
MOS field into four equal quadrants. 

liptical morphology decreases at z > 0.8 (Figure 
[T4| . This decrease is seen with both classifications. 
The decrease is much faster if we consider the C09 
classification which is more conservative in select- 
ing pure elliptical galaxies. This trend could show 
that: the fraction of Spi/Irr with a quenched SFR 
increases with redshift ; the "red and dead" local 
elliptical galaxies didn't have fully acquired their 
elliptical morphology at z > 0.8. This result is 
discussed in Bundy et al. (2009) and Oesch et al. 
(2009). 

6.3. Stellar mass density of quiescent and 
red elliptical galaxies 

The stellar mass density quantifies the total 
stellar mass locked up in a given population (p = 
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Fig. 14. — Fraction (in %) of "quiescent" galaxies 
with an elliptical morphology. The fraction is ob- 
tained by dividing the MF of the red ellipticals by 
the MF of the quiescent galaxies. The upper and 
lower limits are obtained using the G-C and C09 
morphological classifications, respectively. 

4>{M)dM). 

The stellar mass density of quiescent galaxies 
(shown in Figure [15]) increases by 1.1 dex (a factor 
of 14) between z = 1.5 — 2 and z = 0.8 — 1, and 
still increases by 0.2 dex between z = 0.8 — 1 and 
z = 0.2 — 0.4 (by 0.3 dex if we consider the local 
measurement by Bell et al. 2003). Therefore, the 
stellar mass assembly of quiescent galaxies appears 
to slow down at z < 1. 

Figure [15] shows the stellar mass density of red 
elliptical galaxies. The lower and upper limits of 
the shaded areas correspond to the C09 and G- 
C morphological classifications, respectively. At 
z < 0.8, the total stellar mass of red elliptical 
galaxies contributes more than 80% (70%) of the 
total stellar mass in quiescent galaxies. 

6.4. Comparison with the literature 

Figure [TB] compares our "quiescent" MFs with 
those from the literature. The K-band luminos- 
ity functions from Cirasuolo et al. (2007) and 



Fig. 15. — Stellar mass density of quiescent galax- 
ies (black solid circles) and red ellipticals (red open 
circles). The lower and upper envelopes of the 
shaded area correspond to the stellar mass density 
computed using the C09 and G-C classifications, 
respectively. The open star is the local measure- 
ment by Bell et al. (2003). The vertical boxes 
for the quiescent correspond to the two extreme 
values measured in 4 quadrants of 0.5-deg 2 . 

Arnouts et al. (2007) are converted into stellar 
mass functions using the mass-to-light ratio given 
in Arnouts et al. (2007). 

Most of the MFs from the literature are derived 
for red sequence galaxies. Therefore, we split the 
COSMOS sample into "red sequence" and "blue 
cloud" galaxies according to the empirical limit 
M NUV - M R = 0.5 log(M) - 0.8 z - 0.5 (see Fig- 
ure [20]). Our MFs for red sequence galaxies (open 
triangles in Figure fTH]) are systematically higher at 
z > 1 than our MFs for quiescent galaxies (open 
circles). We interpret this difference as being due 
to the presence of dust-extincted spirals within the 
red sequence (see also Figure [5} ■ This contamina- 
tion of the red sequence could increase with red- 
shift since the density of dusty star-forming galax- 
ies increases with redshift (e.g. Takeuchi et al. 
2007, Le Floc'h et al. 2005). 

The agreement between the various surveys at 
the high-mass end is good (within 0.2 dex), given 
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Fig. 16. — Comparison between our MFs derived for quiescent (black open circles) and red sequence (open 
triangles) galaxies and the MFs from the literature. The MFs measured by Arnouts et al. (2007) (solid red 
curves) are based on a template-fitting classification. The MFs measured by Cirasuolo et al. (2007) (green 
dotted lines) are obtained using a red sequence classification. 



the uncertainties in the stellar masses, the dif- 
ferent classification methods and cosmic variance. 
The increase of massive red galaxies between z = 



1.5 — 2 and z — 0.8 — 1 has also been seen in pre- 
vious studies (e.g. Cirasuolo et al. 2007, Arnouts 
et al. 2007). However, the amplitude of this in- 
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Fig. 17. — MF by spectral type. The sample is split into "high activity" (blue vertical shaded area), 
"intermediate activity" (green oblique shaded area) and "quiescent" (red horizontal shaded area) galaxies. 
The lower and upper limits of the shaded areas are the extreme values of the MFs estimated in four quadrants 
of 0.5-deg 2 , which quantifies the impact of cosmic variance. The blue short-dashed lines, the green dotted lines 
and the red long-dashed lines are the MFs measured at z = 0.2 — 0.4 for the "high activity", "intermediate 
activity" and "quiescent" galaxies, respectively. 



crease differs between the samples. At z = 1.5 — 2, 
the MF normalization derived by Cirasuolo et al. 



(2007) is a factor ~10 higher than Arnouts et al. 
(2007). Our MFs for red sequence (open triangles) 
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and quiescent galaxies (open circles) show that the 
two measurements can be partly reconciled when 
consistent selection criteria are applied (Cirasuolo 
et al. used red sequence galaxies and Arnouts et 
al. used a template fitting classification close to 
our quiescent definition). 

Star-forming galaxies 
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Fig. 18. — MF of "star-forming" galaxies (sum 
of intermediate and high activity galaxies) from 
z — 2 to z = 0.2. The vertical box quantifies the 
cosmic variance at z = 0.2 — 0.4 (Scoville et al. 
2007). 

7. Stellar mass assembly of star-forming 
galaxies 

This section presents the MF and stellar mass 
density for star-forming galaxies. We subdivided 
the star-forming sample into "intermediate activ- 
ity" and "high activity" galaxies using a best-fit 
template procedure (see §3.21 and Figure [9]). We 
did not attempt to introduce a morphological sep- 
aration since it would require a too fine classifica- 
tion within the Spi/Irr population. 

7.1. MF of star- forming galaxies 

Figure [17] shows the MF evolution of the "in- 
termediate activity" (green oblique shaded area) 
and "high activity" galaxies (blue vertical shaded 
area). The MF of the "quiescent" galaxies (red 



horizontal shaded area, see §5]) is added as a ref- 
erence. Both the MFs of "intermediate activity" 
and "high activity" galaxies evolve between z = 2 
and z — 0.2. Since their shapes change with time, 
this evolution is mass-dependent. 

As a consistent trend at all redshifts, the slope 
of the "high activity" galaxies is always the steep- 
est. The density of "high activity" galaxies de- 
creases with cosmic time but the size of this de- 
crease is a strong function of the stellar mass. Be- 
tween z = 1.2—1.5 and z = 0.2 — 0.4, their number 
density decreases by a factor of 5 at log(A4) > 11, 
and only by a factor of 1.1 at 9.5 < log(Ai) < 10. 

The "intermediate activity" MF follows a dif- 
ferent evolution. Galaxies as massive as log(M) ~ 
11.6 are already in place at z = 1.5 — 2. The den- 
sity of lower mass galaxies rises with times. 

We plotted in Figure [18] the MFs of all star- 
forming galaxies, i.e. the sum of the "intermediate 
activity" and "high activity" MFs. The exponen- 
tial cutoff does not evolve by more than 0.2 dex be- 
tween z = 1.2 — 1.5 and z = 0.2 — 0.4, which is con- 
sistent with cosmic variance and systematic errors 
in stellar mass measurements. We do not observe 
significant changes in the MF shape. Therefore, 
the decrease with time of "high activity" galaxies 
is partly counter-balanced by the build up of the 
"intermediate activity" MF. 

7.2. Stellar mass density of star- forming 
galaxies 

Figure [T9l shows the evolution of the integrated 
stellar mass density for "intermediate activity" 
(green triangles), "high activity" (blue squares) 
and all star-forming galaxies (sum of intermediate 
and high activity: open circles). 

The stellar mass density increases between z = 
1.5 — 2 and z = 0.8 — 1 for star-forming popula- 
tions. The density increases by 0.5 dex, 0.2 dex 
and 0.3 dex for the "intermediate activity", the 
"high activity" and all star-forming galaxies, re- 
spectively. 

The "intermediate activity" and "high activity" 
galaxies follow an opposite evolution at z < 1 . Be- 
tween z = 0.8—1 and z = 0.2 — 0.4, the stellar mass 
density increases by 0.3 dex for the "intermediate 
activity" galaxies, while it decreases by 0.4 dex 
for the "high activity" galaxies. The stellar mass 
density of all star forming galaxies shows little evo- 
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Fig. 19. — Stellar mass density of star-forming 
galaxies. The green open triangles and the blue 
solid squares are the "intermediate activity" and 
"high activity" galaxies, respectively. The black 
open circles are all star forming galaxies (sum of 
"high activity" and "intermediate activity"). The 
open star is the local measurement by Bell et al. 
(2003). The boxes correspond to the two extreme 
values measured in 4 quadrants of 0.5-deg 2 . 

lution at z < 1. Indeed, the stellar mass density 
at z ~ 0.1 measured by Bell et al. (2003) for star- 
forming galaxies is consistent with the stellar mass 
density we measure at z ~ 1. 

8. Star-forming galaxies with an elliptical 
morphology 

As shown in Figure [20] the elliptical galaxies 
preferentially have a red color. However, a signif- 
icant population of elliptical galaxies appears to 
be blue. The presence of "blue elliptical" is not 
unexpected. The "blue elliptical" galaxies could 
be newly formed elliptical galaxies still consum- 
ing the gas of their progenitors. Accretion of new 
cold gas into an old elliptical galaxy (Hammer et 
al. 2007) could also produce a blue color. We 
present here the MF of this elliptical population 
which is star-forming. 

Figure [21] shows the MFs of the elliptical galax- 



Fig. 20.— Rest-frame color M(NUV) - M{r+) 
(not corrected for dust reddening) as a function 
of stellar mass at i + < 24 in a given redshift bin 
0.9 < z < 1. The red open circles are the galaxies 
morphologically selected as elliptical by C09. The 
green open squares are selected as elliptical with 
the G-C parameters and not by C09. The blue 
crosses are Spi/Irr galaxies classified with the G-C 
parameters. The magenta dashed line is the limit 
adopted to split the sample into red sequence and 
blue cloud galaxies. 

ies purely selected by morphology (black vertical 
shaded area), the "blue elliptical" galaxies (blue 
shaded area) and the red elliptical galaxies (red 
shaded area). The shape of the "blue elliptical" 
MFs differs from those of the red elliptical galax- 
ies: the slope is steeper and the exponential cutoff 
is shifted to lower mass (in agreement with Ilbert 
et al. 2006a). As a consequence of these different 
shapes, the contribution of the "blue elliptical" 
galaxies to the total elliptical population depends 
strongly on stellar mass. Figure [22] shows that, 
regardless of redshift, the fraction of "blue ellipti- 
cals" decreases towards high mass systems. The 
"blue elliptical" galaxies represent < 20% of the 
massive elliptical galaxies (at log(Ai) > 11 and 
z < 1), but their contribution reaches 40-60% at 
log{M) ~ 10. 
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Fig. 21. — MF of elliptical galaxies selected in morphology without any cut in color (black vertical shaded 
area). The blue oblique shaded areas and the horizontal red shaded areas are the MFs of the blue and red 
elliptical galaxies, respectively. The upper and lower envelope MFs are obtained using the G-C and C09 
morphological classifications, respectively. For clarity in the figure, the non-parametric estimates are shown 
only for the C09 classification. The dashed line is the MF of elliptical galaxies derived at z = 0.2 — 0.4, 
which is shown in each rcdshift bin to serve as a reference. 



9. Discussion 

9.1. Redistribution of the star formation 
activity along cosmic time 

We first discuss our results on the MF of star 
forming galaxies. 
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Fig. 22. — Fraction (in %) of blue galaxies in the 
elliptical sample (morphologically selected). The 
upper and lower limits are obtained using the G- 
C and C09 morphological classifications, respec- 
tively. 

The total stellar mass of "high activity" galax- 
ies decreases by a factor of 1.7 between z = 1.5 — 2 
and z = 0.2 — 0.4 (see S}7|), while the stellar mass 
of a given star-forming galaxy grows. A simple 
interpretation is that, after an intense star for- 
mation activity period, these galaxies evolve into 
less active systems (intermediate activity or quies- 
cent). Based on the BC03 models, a passive evo- 
lution could transform a "high activity" galaxy at 
z ~ 1.5 into an "intermediate activity" galaxy at 
z ~ 0.70. 

Figure [521 shows the fraction of "high activ- 
ity" galaxies within the star-forming sample. This 
fraction does not decrease uniformly at all stellar 
masses. Between z = 1.5 — 2 and z = 1 — 1.2, the 
fraction of "high activity" galaxies at log(Ai) ~ 11 
drops by a factor of ~4 (from 40% to 10%), while 
the fraction stays approximately constant at low 
mass log(Ai) ^9.5—10. The contribution of low 
mass, high activity galaxies starts to decrease sig- 



5 A galaxy with a rest-frame color (NUV — r + ) temp [ ate ~ 
1 reddens by 2 mag in 3 Gyr assuming an exponentially 
decreasing SFR with r = 1 Gyr and a solar metallicity. 
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Fig. 23. — Fraction in % of "high activity" galaxies 
within the star-forming sample as a function of the 
stellar mass and per redshift bin (MF of the "high 
activity" divided by the MF of all star-forming 
galaxies). The thick blue line corresponds to the 
stellar mass range covered by the data. The thin 
green line is obtained using the extrapolation of 
the MF with a Schechter function. 

nificantly only at z < 1. Since each spectral type 
corresponds to a range of SSFR (see §3.2|) . it im- 
plies that the low mass star-forming galaxies are 
able to maintain a high SSFR, while the massive 
galaxies evolve rapidly into systems with a lower 
SSFR. Therefore, the redistribution of the star for- 
mation activity follows a clear "downsizing" pat- 
tern (Cowie et al. 1996) within the star-forming 
sample itself. 

The MF of all star-forming galaxies (sum of 
intermediate and high activity) shows little evo- 
lution between z ~ 1.5 and z <~ 0.2. The re- 
distribution of the star formation activity between 
"intermediate activity" and "high activity" galax- 
ies does not affect the overall mass distribution 
of the star- forming galaxies. The little evolution 
of the star-forming MF means that a fraction of 
star-forming galaxies is transferred to the quies- 
cent population (as already noted by Arnouts et 
al. 2007 and Cowie et al. 2008), since star forming 
galaxies generate new stellar populations between 
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z *~ 1.5 and z ~ 0.2. Using backward evolution 
models, Boissier et al. (2009) discuss in detail 
the consistency between the little evolution of the 
star-forming MF and the building of the quiescent 
population. We discuss in the next sections which 
processes can generate these quiescent systems. 

9.2. Mass-dependent assembly of elliptical 
galaxies at z < 1 

We find that the most massive quiescent galax- 
ies are already in place at z ~ 1 while their density 
still rises at low/intermediate masses at z < 1 (see 
This mass-dependent evolution of quiescent 
galaxies confirms the "downsizing" pattern found 
by the COMBO-17 survey (Borch et al. 2006) and 
the DEEP2 survey (Bundy et al. 2006, Cimatti et 
2006). We also investigated the MF evolution of 
the quiescent galaxies with an elliptical morphol- 
ogy. We retrieved the same downsizing pattern: 
the most massive red elliptical galaxies are already 
in place at z ~ 1, while the low/intermediate mass 
red E/S0 galaxies are still being created at z < 1. 
In principle, the weak evolution of the massive red 
elliptical galaxies at z < 1 could be explained by 
a selection procedure based on the galaxy spectra 
(Van Dokkum & Franx 2001) since blue elliptical 
galaxies missed by a multi-color criterion could ac- 
count for a significant evolution of the high-mass 
end. However, we showed in section that the 
contribution of blue elliptical galaxies is limited at 
20% at log(M) > 11. Therefore, the blue elliptical 
contribution can not explain the low evolution rate 
of the most massive elliptical galaxies at z < 1. 

A possible interpretation of this mass depen- 
dent evolution is a galaxy assembly by merger pro- 
cess more efficient at low/intermediate mass than 
at high mass at z < 1, if we assume that red ellip- 
tical galaxies are formed by mergers (e.g. Toomre 
& Toomre 1972, Athanassoula 2008, Bekki et al. 
2008)@. This picture is in agreement with De Ravel 
et al. (2008) who found that the merger rate de- 
creases with stellar mass when using galaxy pair 
counts in the VIMOS-VLT Deep Survey (le Fevre 
et al. 2005). 

The formation of low mass quiescent galaxies 



Bckki ct al. (2008) simulated merging between gas rich 
spiral galaxies of mass JV[ ~ 10 9 .Mq, which formed a dy- 
namically relaxed low mass elliptical galaxy in 1.4 Gyr with 
a SFR of 0.03 M Q /yr. 



with an elliptical morphology could also be ex- 
plained by "morphological quenching" (Martig et 
al. 2009). The presence of a spheroid is suffi- 
cient to stabilize the gas disk and quench the star 
formation. This process is efficient in low mass 
halos which could explain the formation of these 
low mass, quiescent and elliptical galaxies. 
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Fig. 24. — Evolution of the stellar mass density 
of star-forming (blue filled circles) and quiescent 
galaxies (red filled circles) compared to various re- 
sults from the literature. The stellar mass density 
of all the red sequence and quiescent galaxies is 
shifted vertically by -0.5 dex for the clarity of the 
figure. The offset of 0.2 dex between Arnouts et 
al. (2007) and our measurement for star-forming 
galaxies is discussed in appendix [Pi 

9.3. Slow down in the assembly of massive 
quiescent galaxies at z ~ 1 

The most massive quiescent galaxies are not in 
place at z ~ 2 and their number density evolves 
rapidly between z ~ 2 and z ~ 1. The exponen- 
tial cutoff of their MF increases by 0.4 dex between 
z = 1.5 - 2 and z = 0.8 - 1. Figure |2H shows that 
the rapid assembly rate of quiescent galaxies at 
1 < z < 2 is consistent with other surveys, de- 
spite the different methods used to classify early 
type galaxies, to determine the stellar masses and 
to measure the distances (spectro-z or photo-z). 
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Fig. 25. — Fraction of the "quiescent" over the 
total population per stellar mass bin (MF of the 
quiescent divided by the total MF). The thick red 
line corresponds to the stellar mass range covered 
by the data. The thin green line is obtained using 
the extrapolation of the MF with the Schechter 
function. 

Using deeper data in the GOODS field, we also 
checked that our results don't suffer from a signif- 
icant incompleteness which could mimic the rapid 
assembly of quiescent galaxies at 1 < z < 2 (see 
appendix |EJ . 

Therefore, the most massive quiescent galaxies 
are created rapidly between z ~ 2 and z ~ 1 and 
their assembly slows down at z < 1. We tenta- 
tively explain this slow down in their evolution by 
analyzing the relative evolution of the quiescent 
and star-forming MFs from z = 2 to z = 0.2. 

Figure [231 shows the fraction of quiescent galax- 
ies as a function of stellar mass. The quies- 
cent galaxies represent less than 20% of the most 
massive galaxies at z — 1.5 — 2. As a con- 
sequence, "wet mergers" between massive star- 
forming galaxies directly create new massive quies- 
cent galaxies (since star-forming galaxies are more 
massive and more numerous) which generates a 
rapid growth of the quiescent high mass end be- 
tween z ~ 2 and z ~ 1. 

A change of regime occurs at z < 1, where 



the exponential cut-off of star-forming galaxies is 
shifted at lower mass than the exponential cut-off 
of quiescent galaxies (see Figure [TT]) . This change 
is likely explained by the rapid increase of the mas- 
sive quiescent population at 1 < z < 2, combined 
with the decrease of the SFR at z < 1 (e.g. Lilly 
et al. 1996, Le Floc'h et al. 2005, Tresse et al. 
2007) which prevents to regenerate the massive 
star-forming population. Therefore, the quiescent 
population dominates the massive end of the MF 
at z < 1, with more than 60-70% of quiescent 
galaxies at log(A4) > 11 — 11.5. As a consequence, 
"wet mergers" become inefficient at z < 1 to gen- 
erate the most massive quiescent galaxies. 

Dry merging (merging between quiescent galax- 
ies, e.g. Van Dokkum et al. 2001) is the only 
process left to form the most massive quiescent 
galaxies at z < 1. This process has less impact on 
the quiescent MF evolution since "dry merging" is 
not a direct supply of new quiescent systems (the 
progenitors are already quiescent galaxies). More- 
over, "dry merging" involving two massive quies- 
cent galaxies (log(A4) > 11) is not a common pro- 
cess since the exponential cutoff of the quiescent 
MF does not evolve significantly at z < 1. 

Therefore, the disappearance of "wet merging" 
as an efficient process to form the most massive 
ellipticals could explain the slow down in their as- 
sembly at z < 1. 

9.4. Indirect constraints on the AGN feed- 
back 

Two different modes of AGN activity are usu- 
ally considered: the "bright mode" (e.g. QSO) 
and the "radio mode" which radiates less energy 
than the "bright mode" but is a more common 
mechanism (e.g. Croton et al. 2006). In both 
modes, the AGN feedback can prevent the gas 
to cool in the dark matter halo. Since the "ra- 
dio mode" is less energetic, the radio jet from the 
black hole prevent the gas to cool only in quasi- 
hydrostatic shock-heated halos, i.e. in halos above 
a critical mass around 10 12 M Q (Cattaneo et al. 
2006, Somerville et al. 2008). 

As a consequence, the radio mode prevents the 
star formation of recurring in all the galaxies of 
massive halos. Therefore, the star formation is 
shut down in the galaxy hosting the AGN, but 
also in all the galaxies surrounding the AGN. We 
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Fig. 26. — MF of "quiescent" galaxies (thick red 
lines in the top panel) and elliptical galaxies (thick 
red lines in the bottom panel) at z = and z — 1 
predicted by Cattaneo et al. (2008) using the 
GallCS semi-analytical simulations and a "halo 
quenching" model. Our quiescent and red ellip- 
tical MFs are reported in each panel. 

discuss here the impact of AGN feedback on the 
"surrounding galaxies"0- A specific study on the 
host galaxies of the radio sources is presented in 
Smolcic et al. (2009). 

The shutdown of the star formation in the "sur- 
rounding galaxies" does not modify their morphol- 
ogy. As their original morphology can be any, we 
then expect a significant population of quiescent 
galaxies with a Spi/Irr morphology. Indeed, the 
ratio between the red elliptical and the quiescent 
MFs shows a significant fraction of Spi/Irr hav- 
ing a quenched star formation at low masses (40- 
60% at log(M) ~ 9.5, see Figure[l4|), which leaves 
room for this "external" action mode of AGN feed- 
back. 

More quantitatively, Cattaneo et al. (2008) 
provided the predicted MFs for quiescent (red se- 



We call "surrounding galaxies" the galaxies which belong 
to the same halo as the AGN but which are not hosting 
the AGN themselves. These galaxies are satellite galaxies 
if the AGN is hosted by the central galaxy of the dark 
matter halo. 



quence in Cattaneo et al.) and elliptical (bulge- 
dominated in Cattaneo et al. ) galaxies at z = 
and = 1. They used a "halo quenching" model 
(Somerville et al. 2008). In this model, the 
AGN feedback shuts down the star formation in 
all the galaxies within a halo more massive than 
10 12 M Q (Cattaneo et al. 2006, Somerville et al. 
2008). We report in Figure [26] the MF predicted 
by Cattaneo et al. (2008). As expected, the pre- 
dicted density of quiescent galaxies stays constant 
at low/intermediate masses (top panel), while the 
density of elliptical galaxies decreases toward low 
masses (bottom panel). On the basis of these 
predicted MFs, ~ 10% of the quiescent galaxies 
would be elliptical at log(M) ~ 10. By contrast, 
we found a fraction of 80%(60%) quiescent galax- 
ies with an elliptical morphology (see Figure [14]) . 
Therefore, quenched star formation is more often 
linked to an elliptical morphology than would be 
predicted by a "halo quenching" model. These 
measurements constrain in the same way all the 
processes quenching star formation without mod- 
ifying the galaxy morphology, for example, gas 
starvation (i.e. satellite galaxies not fueled in cold 
gas since they are not at the center of the dark 
matter halo potential well). 

A more quantitative constraint on these mech- 
anisms requires a detailed comparison with pre- 
dictions of semi-analytical simulations, which is 
planned for forthcoming papers. 

10. Conclusion 

We derived the galaxy stellar mass function and 
stellar mass density in the 2-deg 2 COSMOS field. 
We explored stellar mass assembly by morpholog- 
ical and spectral type from z = 2 to z — 0.2. The 
MF estimate is based on 196,000 galaxies selected 
at F3.6^m > lfJ-Jy and photo-z with an accuracy 
of CT ( Zp w-^c)/(i+* 3P ec) = 0-008 at i+ B > 22.5. 
We summarize our results below. 

• We found that z ~ 1 is an epoch of transi- 
tion in the assembly of quiescent galaxies. Their 
stellar mass density increases by 1.1 dex between 
z = 1.5 — 2 and z — 0.8 — 1 (corresponding to a 
period of 2.5 Gyr), but only by 0.3 dex between 
z = 0.8—1 and z ~ 0.1 (a period of 6 Gyr). 
The high- mass exponential cutoff of the quiescent 
MF increases by 0.4 dex between z = 1.5 — 2 and 
z = 0.8 — 1, but almost no evolution is seen at 
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z < 1. We investigated if the weak evolution of the 
most massive quiescent galaxies is also seen using a 
morphological classification. The exponential cut- 
off of massive red elliptical galaxies does not in- 
crease significantly at z < 1. Moreover, the blue 
elliptical galaxies do not contribute more than 20% 
to the high-mass end of the total elliptical sample, 
which is not sufficient to produce significant evo- 
lution in the exponential cutoff. 

• We found that the high-mass end of the star- 
forming MF is shifted below the high-mass end of 
the quiescent MF at z < 1. Therefore, we inter- 
preted the slow down in the assembly of the most 
massive elliptical galaxies at z < 1 as being due to 
a "lack of supply" of massive star-forming galaxies 
available for "wet mergers" . 

• We observed a rapid rise of quiescent galaxies 
at low/intermediate masses. We characterized the 
nature of these newly formed quiescent galaxies by 
adding morphological information. We quantified 
the fraction of quiescent galaxies with an ellipti- 
cal morphology, as well as the fraction of Spi/Irr 
galaxies with quenched star formation. The sig- 
nificant fraction of quenched Spi/Irr (40 — 60% at 
log(M) ~ 9.5) leaves room for a mechanism which 
shuts down the star formation without transform- 
ing their morphology, such as the impact of AGN 
feedback on the satellite galaxies of a massive halo 
(e.g. Cattaneo et al. 2006). Since the majority 
of quiescent galaxies have an elliptical morphol- 
ogy at z < 0.8 (80-90% at log(M) ~ 11), the 
dominant process which shuts down star formation 
should be linked to the acquisition of an elliptical 
morphology, as might be expected in galaxy merg- 
ing and/or morphological quenching (Martig et al. 
2009). 

• Finally, we divided the star-forming sample 
into "intermediate activity" and "high activity" 
galaxies, which corresponds to two classes of SSFR 
(SFR divided by stellar mass). The MF of the 
"high activity" galaxies shows that the most mas- 
sive of them end their high activity phase first. 
Therefore, the low mass star-forming galaxies are 
able to maintain a high SSFR, while the massive 
galaxies evolve rapidly into systems with a lower 
SSFR. This redistribution of the star formation ac- 
tivity follows a clear "downsizing" pattern (Cowie 
et al. 1996) within the star-forming sample itself. 
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Table 2: Schechter parameters of the MFs for the quiescent and elliptical galaxies between z = 0.2 and z = 2. 
The errors combined the one sigma Poissonian errors (2AlnC — 1) as well as the uncertainties induced by 
the photo-z. Parameters a listed without errors are set "ad-hoc", and the errors on log(M*) and are 
obtained by varying a by ±0.5 around the fixed value. The morphological classification is obtained using 
the G-C parameters. However, we caution the reader that the errors are probably underestimated: the error 
budget is dominated by systematic effects (e.g. possible photo-z biases, systematic uncertainties related to 
the stellar mass estimate) that are not included here. 
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Table 3: Schechter parameters of star-forming galaxies ("intermediate activity", "high activity") between 
z = 0.2 and z = 2. A parametrization of the total MF can be retrieved by summing the Schechter fit of the 
"quiescent" (given in Table [2]), "intermediate activity" and "high activity" galaxies. Errors are computed as 
in Table 2. 
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A. Absolute magnitude estimate 

At high rcdshift, the k-correction is one of the main sources of systematic error in the absolute magnitude 
and rest-frame color estimate. The k-correction depends on the galaxy SED, which is not directly observed 
(Oke et Sandage 1968). In order to minimize the uncertainty induced by the k-correction term, the rest- 
frame luminosity at a given wavelength, A, is derived from the apparent magnitude observed at A x (1 + z) 
(appendix A of Ilbert et al. 2005). With this procedure, the absolute magnitudes are less dependent on the 
SED. One drawback of this method is that the uncertainty in the observed apparent magnitude is directly 
propagated into the absolute magnitude. For this reason, the absolute magnitudes were measured from one 
of the following bands u*g + r + i + z + K and 3.6/im all of which have the highest signal-to-noise ratio and a 
zero-point correction lower than 0.05 mag (Ilbert et al., 2009). 



B. Additional constraint in the template-fitting using the 24/ito fluxes 

The deep MIPS S-COSMOS data were taken during Spitzer Cycle 3 and cover the full COSMOS 2-deg 2 
(Aussel et al. 2009, in prep.). The 24/ito sources were detected with SExtractor (Bertin & Arnouts 1996) 
and their fluxes measured with a PSF fitting technique (Le Floc'h et al. 2009). The infrared luminosities 
L^ IPS were extrapolated from the 24 \xm fluxes using the Dale & Helou (2002) library (Le Floc'h et al. 
2009) and converted into SFRin using the calibration from Kennicutt (1998). 

The stellar mass estimate and instantaneous SFR (hereafter S F R tem piate) are derived from the best-fit 
template. We compared SFR tem piate and SFRir in order to decide which extinction law to adopt - either 
Calzetti et al. (2000) or Chariot & Fall (2000). Figure [27] shows that the use of the Calzetti et al. (2000) 
extinction law reduces the systematic offset between S F Rtempiate and SFRir. Therefore, we favored the 
Calzetti et al. extinction law for our analysis. 

As a second step, we used the 24 /im MIPS fluxes as an additional constraint in the template fitting pro- 
cedure. The goal of this additional constraint is to remove possible degeneracies between "old and quiescent" 
models and "star- forming and dust extincted" models. For each template, we computed i}^™v late ; which is 
the infrared luminosity which would be re-emitted in the infrared according to the template, assuming that 
all of the UV light absorbed by dust is re-emitted in the infrared and that the massive stars are the only 
source of infrared emission. A likelihood, £, is computed at each step of the template fitting procedure. We 
multiplied this likelihood by the probability to measure L™ IPS for a given template: 

/tMIPS t template\2 

CXp( R ~ 1B 2 — ) 

C' oc £(template, scaling, E(B — V)) X ^II_ (Bl) 

err x y2n 

where err is the error on the L^ IPS measurement (Le Floc'h et al. 2009, in prep.). We added 0.2 dex 
in quadrature to the error on Lf^ IPS to take into account systematic uncertainties in L^ IPS . When the 
galaxy has no MIPS counterpart at more than 2", an upper limit is applied to the °* e _ This upper- 

limit corresponds to the lowest L^ PS observable for a sample selected at Fi^ m > 100/iJy (5er completeness 
limit). No prior is applied if the optical counterpart is between 0.6" and 2" (less secure optical counterpart). 

By construction, the comparison between SFRtempiate and SFRir is improved by this prior as shown in 
Figure EH 



C. The tool ALF 

The selection of the galaxy sample at Fs.e^m > l^Jy defines a limit in rcdshift above which the galaxies 
are too faint to be observed. Statistical estimators are required to correct for the incompleteness created by 
this flux limit. 

We measured the luminosity and stellar mass functions using the tool ALF (Algorithm for Luminosity 
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0.5 1 1.5 2 2.5 

SFR from mid-IR (24jum) 



Fig. 27. — SFR estimated from the best-fit template versus the SFR directly measured using the mid-infrared 
24 /j,m MIPS flux. The short-dashed blue contours and solid black contours are obtained using the Chariot 
& Fall (2000) extinction law and the Calzetti et al. (2000) extinction laws, respectively. The long-dashed 
red contours are obtained using the constraint on the 2A(im flux. 

Function) described in Ilbert et al. (2005) and originally developed to measure the luminosity functions from 
the VIMOS-VLT Deep Survey (le Fevre et al. 2005). The procedures used to compute the stellar MF or the 
luminosity function are the same. The tool includes various estimators: the non-parametric 1/V maxi C' + , 
SWML and the parametric STY. The STY (Sandage et al. 1979) and SWML (Efstathiou 1988) determine 
the MF by maximizing the likelihood to observe a given stellar mass-redshift sample. The STY estimator 
presupposes that the MF can be parametrized with a Schechter function (Schechter 1976): 

$>(M)dM = <f>*(M/M*) a exp(-M/M*) d(M/M*). (CI) 

This parametrization allows us to describe the MF using three parameters: a (slope), Ai* (characteristic 
stellar mass) and </>* (normalization). The SWML is a non-parametric estimate of the MF, useful to verify 
that a Schechter function is a good representation of the data. The non-parametric 1/V max estimator 
(Schmidt 1968) is the most widely used because of its simplicity. The 1/V max is the inverse sum of the 
volume in which each galaxy could be observed. The 1/V max is the only estimator directly normalized. 
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Lynden-Bell (1971) derived the non-parametric C~ method to overcome the assumption of a uniform galaxy 
distribution derived using 1/V max (we used a slightly modified version called C + , Zucca et 1997). The 
implementation of these estimators is detailed in appendix A of Ilbert et al. (2005). 



9.0 < log(M) < 9.3 10.5 < log(M) < 11.5 




0.5 1 1.5 2 -0.5 0.5 

photo — Z log(stellar mass/L^) - Arnouts et al. 2008 



Fig. 28. — Left panel: mass-to-light ratio as a function of redshift. The stellar masses are computed from a 
SED fitting procedure. The luminosity is measured in the rest-frame if s -band. The top and bottom panels 
correspond to the blue cloud and red sequence galaxies, respectively. The solid lines correspond to the 
empirical relations derived by Arnouts et al. (2007). The relation of Arnouts et al. (2007) is derived for the 
massive galaxies brighter than m^^^ m < 21.5. The magenta crosses and the green circles are the galaxies 
selected in stellar mass ranges 9.0 < log(M.) < 9.5 and 10.5 < log(Ai) < 11.5, respectively. Right panel: 
difference between the observed mass-to-light ratio and the analytical relation established by Arnouts et al. 
(2007). The top and bottom panels correspond to the star- forming and quiescent galaxies, respectively. The 
difference is shown per bin of stellar mass. 



D. Mass-to-light ratio 

The if-band luminosity is not a direct tracer of the stellar mass. Numerous analyzes have used an 
analytical parametrization of the mass-to-light ratio in order to derive the stellar masses (e.g. Kochanek et 
al. 2001, Driver et al. 2007, Arnouts et al. 2007). 

The left panel of Figure [28] shows the mass-to-light ratio as a function of redshift (the stellar masses are 
derived as described in section l4~T|) . We observed a decrease of the mass-to-light ratio with redshift, in good 
agreement with the empirical relations derived by Arnouts et al. (2007). We also found that the mass-to- 
light ratio increases with the stellar mass. This trend is expected since the high mass galaxies are older and 
have a larger mass-to-light ratio. The right panel of Figure [28] shows the difference between the observed 
mass-to-light ratio and the analytical relation derived by Arnouts et al. (2007). For the quiescent galaxies, 
the median difference is below 0.1 dex. The median difference is also below 0.1 dex for the star-forming 
galaxies (log(M) > 10.5), but the median difference is shifted below -0.3 dex at low masses log(A4) < 9.5, 
showing that this relation is not anymore a good proxy for the low mass galaxies. 

We compared our stellar mass densities for star-forming galaxies with results from the literature in Fig- 
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ure [24] The stellar mass densities of Arnouts et al. (2007) are systematically 0.2 dex higher than our 
measurements. Arnouts et al. (2007) based their measurements on the iC-band LF converted into stellar 
mass density using the mass-to-light relation discussed above. At low masses (log(A4) < 10), we found 
mass-to-light ratio lower by 0.2-0.4 dex than those estimated from the Arnouts et al. (2007) relation. Since 
star-forming galaxies are mostly low/intermediate mass galaxies (steep slope of their MF), the dependency 
of the mass-to-light ratio on the mass likely explains this offset of 0.2 dex in the stellar mass densities. 

E. The quiescent MFs in the GOODS and COSMOS fields 

Two possible sources of incompleteness could affect our estimate of the quiescent MFs at z > 1: the depth 
of the optical catalogue and the confusion in the IRAC images. We checked the validity of our selection 
criteria using the GOODS datasets (Giavalisco et al. 2004). We used the public catalogues FIREWORKS 
(Wuyts et al. 2008) and MUSIC (Santini et al. 2009). The IRAC photometry in these catalogues is less 
affected by the confusion since specific softwares allow an accurate "PSF-matching" (e.g. convphot in Santini 
et al. 2009). The GOODS optical data are also deeper than COSMOS (90% complete in ACS/F775W at 
26.5 mag). 

We selected the FIREWORKS and MUSIC catalogues at F(3.6/im) > IfiJy, which is the selection limit of 
our study. We computed the photo-z using the code "Le Phare" and the same setting (templates, extinction, 
...) as Ilbert et al. (2009). The rest-frame colors and stellar masses were estimated following exactly the 
method described in section 14.31 We also classified galaxies according to their dust-corrected NUV — R 
rest-frame colors as in section I3~2"1 

First, we find that the contribution of galaxies fainter than / = 26.5 is negligible at the considered stellar 
mass limits (see tabled]). Then, we derive the quiescent mass functions in GOODS (using both MUSIC and 
FIREWORKS catalogues, as well as their original photo-z). The comparison with the S-COSMOS MFs is 
shown in Figure (29] The results are in excellent agreement. The differences could be easily explained by 
cosmic variance since the GOODS field covers only 160 arcmin 2 . The only significant difference is seen using 
the GOODS-MUSIC catalogue and their own set of photo-z at z > 1.2 (a factor ~ 1.5 in normalisation). 
From this comparison, we conclude that our results don't suffer from any significant incompleteness at z > 1 
and that the rapid assembly of the quiescent population at 1 < z < 2 can't be explained by incompleteness 
or confusion effects. 
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Fig. 29. — Stellar mass functions of the quiescent population in the GOODS and COSMOS fields. The black 
filled circles correspond to the S-COSMOS estimate of the quiescent MFs. The open points are the quiescent 
MFs measured in the GOODS field. Red and green circles are an estimate based on the FIREWORKS 
catalogue using "Le Phare" and the Wuyts et al. (2008) photo-z, respectively. Blue and orange triangles 
are an estimate based on the MUSIC catalogue using "Le Phare" and the Santini et al. (2009) photo-z, 
respectively. 
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